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We present the first calculation of the basic properties of the /-mode instability in rapidly rotat- 
ing relativistic neutron stars, adopting the Cowling approximation. By accounting for dissipation in 
neutron star matter, i.e. shear/bulk viscosity and superfluid mutual friction, we calculate the asso- 
ciated instability window. For our specific stellar model, a relativistic polytrope, we obtain a mini- 
mum gravitational growth timescale (for the dominant £ = m = i mode) of the order of 10^ — 10^ s 
near the Kepler frequency while the instability is active above ~ 0.92 and for temperatures 
~ (10^ — 2 X 10^'^) K, characteristic of newborn neutron stars. 



Introduction — A classic result in gravitational physics 
is that pulsation modes in rotating compact stars can 
be driven unstable under the emission of gravita- 
tional radiation. This is accomplished by the so-called 
Chandrasekhar-Friedman-Schutz (CFS) mechanism |1, 
|2] which is based on the general notion of oscillations 
changing from counter- to co-rotating (with respect to 
the stellar rotation) as a result of rotational dragging. 

It is envisaged that CFS-unstable modes could play 
a key role in neutron star astrophysics by producing de- 
tectable amounts of gravitational radiation which can be 
used to probe the neutron star interior structure ||3j|4|. 
Additionally, they could provide a mechanism for lim- 
iting the stellar spin frequency |5|. The realization of 
such gravitational-wave 'asteroseismology' is expected 
to take place with the planned Einstein Telescope or 
more optimistically even with the Advanced LIGO and 
VIRGO instruments in the near future. 

The initial work on the subject of gravitational 
radiation-driven instabilities was focused on the insta- 
bility of the /-mode in Newtonian neutron star models 
||6H2l and concluded that such systems might be observ- 
able sources of gravitational waves, even beyond our 
Galaxy, in the form of newborn neutron stars spinning 
near the mass-shedding limit. Subsequent work in the 
field was mostly driven by the discovery of a similar 
instability in the inertial r-modes, and its potential rel- 
evance for systems like the accreting neutron stars in 
LMXBs |51. 

In more recent years, progress in numerical hydro- 
dynamics has allowed the modelling of rapidly rotat- 
ing relativistic neutron stars. Several calculations ||T0HT4| 
have shown that relativistic stars enter the /-mode in- 
stability region at a lower spin rate as compared to their 
Newtonian counterparts. This is an indication of an en- 
hanced /-mode instability in relativistic stars, and pro- 
vides strong motivation for a renewed investigation of 
the subject. 

With the purpose of corroborating the case of a 
stronger relativistic /-mode instability, in this Letter we 



present the first calculation of the growth rate of this in- 
stability in relativistic stars and the associated 'instabil- 
ity window', that is, the stellar parameter space where 
the mode growth due to gravitational radiation over- 
comes the dissipative action of neutron star matter. 

The f-mode in rapidly rotating relativistic stars — We 
model a neutron star as a self-gravitating relativistic 
fluid obeying a polytropic equation of state. The pulsa- 
tions of this system are studied by a linear time-domain 
numerical code which evolves the general relativistic 
hydrodynamical equations V^iT''" = in the fixed 
spacetime of the unperturbed 'background' star. A key 
property of our model is its ability to treat rapidly ro- 
tating stars, with the stellar angular frequency O only 
restricted by the Kepler mass-shedding limit Ok- This 
numerical scheme has been recently applied to the gen- 
eral study of the oscillation spectrum of rapidly rotating 
(both uniformly and differentially) relativistic stars IITSl - 

More relevant to the scope of this Letter is our recent 
work IfTSl which has provided state-of-the-art results on 
the /-mode oscillation frequency and gravitational radi- 
ation damping/ growth rate as a function of rotation for 
a collection of polytropic models. Referring the reader to 
that paper for details on the numerical techniques, here 
we simply outline the general strategy followed. A con- 
sistent calculation of the mode's gravitational radiation- 
induced damping /growth rate would require to move 
beyond the Cowling approximation and consider the 
evolving spacetime itself. Since this requires very large 
meshes and long evolution times, and is therefore im- 
practical at present, we content ourselves with an ap- 
proximate solution to the problem, using the standard 
multipole expansion for the gravitational wave flux. 
This is (e.g. EI), 

Egw = Tr E ^f'^'^N, (\3Dem\^ + \Skm\^) (1) 

where cry and cr, = (7,. — mVL is the mode frequency 
measured in the rotating and inertial frame respectively. 
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SDg,„ and SJ(y„ are the mode's mass and current multi- 
pole moments, and Nf is the usual coupling constant. 

In calculating the gravitational wave flux we use the 
relativistic mode eigenfunctions in the volume integral 
for SD(„, and neglect the current multipoles here. We 
follow a similar procedure for calculating the mode's to- 
tal energy Ej^ (measured in the rotating frame). Since 
we apply the Newtonian expressions here, though 
weighted using the proper spatial volume element on 
the spacelike hypersurface, we introduce a slight incon- 
sistency into the calculation. The induced error should 
be comparable to using the quadrupole formula for 
gravitational wave extraction |19|, and is expected to be 
small compared to other modeling uncertainties. 

The timescale associated with gravitational wave 
damping /growth is given by (see, for example, [7]) 
fgw = 2Ein/Egw The instability is indicated by a neg- 
ative fgw < 0. In turn, this requires a cr,- < which 
is equivalent to the statement that (1 should exceed in 
magnitude the mode's "pattern speed" — 17>/ m. This is 
a short description of the CPS instability mechanism jJI- 

Besides the action of gravitational radiation we also 
need to account for dissipation. In a neutron star there 
are several viscosity mechanisms that could potentially 
limit or suppress an instability Here we assume a 
mixture of degenerate neutrons, protons and electrons 
where dissipation is due to the familiar shear and bulk 
viscosities. This is strictly true for normal, i.e. non- 
superfluid matter. At temperatures T < 10^ K the likely 
presence of neutron and proton superfluid condensates 
leads to a multi-fluid system which features several addi- 
tional shear /bulk viscosities [20| and the so-called mu- 
tual friction. The latter effect is included in our analysis 
but we will ignore the extra viscous degrees of freedom 
since they are expected to be comparable to the ones of 
the single-fluid system. In this case mutual friction is 
likely the dominant damping mechanism in the super- 
fluid regime. 

The nature of the standard viscosity is well known. 
Shear viscosity is dominated by electron (neutron) colli- 
sions when neutrons are superfluid (normal) ||2T1 . The 
/5-equilibrium reactions responsible for bulk viscosity 
are assumed to be those of the modified Urea process 
Il22l . The associated viscosity coefficients can be found 
in the literature Il7l l2l]l : the same coefficients are used in 
the present work (also, a more elaborate study on shear 
viscosity including the dependency on the proton frac- 
tion does not change our results significantly ESl ). 

The more complex issue of the formation of a viscous 
boundary layer at the crust-core interface has never been 
discussed in the context of the /-mode, although it is 
known to be a key damping mechanism for the r-mode 
instability (see e.g. |5J and references therein). It is ne- 



glected here, firstly, because we are mostly interested in 
newborn neutron stars with temperatures comparable 
to the melting temperature of the crust 10^^ K [22]); 
at this stage the crust physics is highly uncertain and, 
secondly, because in colder stars the instability is most 
likely suppressed by mutual friction (see discussion be- 
low). 

Mutual friction in neutron stars is a less familiar con- 
cept. When the bulk of the interior matter is in a su- 
perfluid state, the possibility of large scale relative flow 
between the fluids is counteracted by mutual friction, 
which refers to the vortex-mediated coupling between 
the fluids. It can be quantified in terms of the drag pa- 
rameter TZ which is (roughly) the ratio of the drag force 
to the Magnus force experienced by a moving neutron 
vortex. The /-mode in superfluid neutron stars generi- 
cally contains a 'counter-moving' degree of freedom and 
is therefore susceptible to damping by mutual friction 

L8,24J. 

In accordance with the strategy followed for the 
growth timescale fgw, we calculate the shear and bulk 
viscosity damping timescales (fgv and f^v/ respectively) 
by using the inviscid relativistic mode eigenfunctions in 
Newtonian formulae (these can be found in fTl)- 

Our single-fluid numerical scheme is oblivious to su- 
perfluid hydrodynamics. Yet, our analysis should some- 
how account for mutual friction since it is likely the 
dominant dissipative effect in 'cold' neutron stars. In 
order to do so, we rely on results from calculations in 
Newtonian superfluid models. From the recent work 
in f24] we can extract the mutual friction damping 
timescale, 

^ (2^ + 3)3(2£ + 5) /c^\ 

fmf~ 2(£ + i)(3^ + 5) \gm) \ nn J ^' 

where Cs is the sound speed and M,R are the stellar 
mass and radius. This expression should be adequate 
for our purposes. 

The combined total dissipation leads to the following 
expression for the mode energy (e.g. f5"7]), 

Em = -2E^ (-^ + (3) 
\ 'gw tdiss / 

where the various dissipation timescales add up as "par- 
allel resistors", f^^ = f-^ + f^J + rj. 

The f-mode instability window — In the absence of dis- 
sipation the /-mode instability is active for any rota- 
tion n > (Ic where the critical rotation is defined by 
(J, (He) = 0. Relativistic stars generally enjoy lower 
critical frequencies (and therefore are more instability- 
prone) as compared to Newtonian models |12J. For the 
relativistic star considered in the following discussion. 
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FIG. 1. The relativistic f-mode instability window. The curves 
are the solutions of — fgw = ^diss' for a N = 0.73 polytropic 
model with M = 1.48 M©, R = 10.47km (in the n = limit). 
The shaded area indicates the presence of superfluidity (for 
a fiducial Ten = 10^ K). The dashed curves represent different 
choices for the mutual friction drag TZ (shown only for the ni = 
4 multipole). 

and for the particular m = A /-mode multipole (which 
is the dominant one as we shall shortly see) the thresh- 
old for the CFS instability is at 0^/0^ w 0.83. 

In the full dissipative system the boundary of the /- 
mode instability region is determined by £m = which 
according to l|3) is equivalent to — fgw = ^diss- The result- 
ing n(T) instability curve is shown in Fig. [T|and is the 
main result of this Letter. The particular example shown 
in the figure corresponds to a neutron star with repre- 
sentative values for mass and radius, approximated by 
aN = 0.73 polytropic model which at the non-rotating 
limit leads to M = 1.48 M© and R = 10.47 km. For this 
model the Kepler frequency is /k = 0x7271 = 1.093 
kHz. 

According to Fig. [T| the m = 4 multipole is the one 
associated with the widest instability window. This 
mode's inertial frame frequency is /, = (Ti/Itz = 
640 — 2290 Hz, corresponding to the instability region 
O/Ok ~ 0.92 — 1. The next lower multipole m = 3 is 
unstable for O/Ok ~ 0.96 — 1, spanning the frequency 
range /, = 200 — 1140 Hz. In contrast, the instability 
window for the m = 2 mode is negligibly small, as it 
is the case in Newtonian calculations [7J. It should be 
pointed out, however, that a finite-size m = 2 instability 
window is present in models with stiffer, but not unre- 
alistic, equations of state. 

As anticipated, the /-mode instability is blocked by 
bulk viscosity in the high temperature regime T > 
2 X 10^*^ K. At lower temperatures the instability has to 
compete against shear viscosity and vortex mutual fric- 
tion. The former effect suppresses the instability below 




FIG. 2. The /-mode growth timescale for the lowest three 
multipoles as a function of the stellar frequency. The stellar 
model is the same as in Fig. [T] The shaded area represents the 
range of the cooling timescale jif , from an initial temperature 
T = 5 X lO^^KtoafinalT = Ten = (5 - 9) x lO^K. 

~ 10'' K. The impact of mutual friction is rather more 
severe. The instability is completely suppressed below 
a critical temperature Ten ^ 10^ K for neutron superflu- 
idity in the core. This result assumes 'standard' mu- 
tual friction, i.e. scattering of electrons by the magne- 
tised neutron vortices, in which case TZ ~ 10~^ Il25l . It 
is possible, however, that mutual friction is dominated 
by the interaction between vortices and quantised flux- 
tubes. This interaction could lead to vortex pinning; for- 
mally, this would correspond to 7^ ^ 1, where it is un- 
clear what the true strength of TZ is. It is likely that the 
assumed mutual friction force leading to l|2} is not appli- 
cable in the case of direct vortex-fluxtube interactions. If 
we treat 7^ as a free parameter, we find that mutual fric- 
tion has a negligible effect on the instability for TZ > 10^. 

Not surprisingly, the overall shape of the relativistic 
instability window does not differ much from the one of 
Newtonian models fT', although it is clearly larger 
What is different is the instability's growth timescale 
fgw Numerical results for the growth rate of the first 
few multipoles are displayed in Fig. |2] Focussing on 
the m = A and m = 3 modes we find a timescale 
fgw ~ 10'^ — 10^ s at n w (Ik- Both modes grow (al- 
most) equally fast above 0.98 (Ik, however the m = 3 
mode may be more detectable due to its lower /. The 
growth timescale is a factor 10 — 100 shorter than the cor- 
responding Newtonian results |4, 7J. This result, which 
can be attributed to the lower relativistic n^, is clearly 
good news for the gravitational wave observability of 
the /-mode. 

The astrophysical systems in which the /-mode in- 
stability is more likely to operate are newborn neutron 
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stars, spinning with Ci > 0.9 Ci^ and with a temperature 
T < 2 X 10^*^ K. Once the star is inside the instability 
region, the efficiency of the /-mode as a source of grav- 
itational radiation is largely decided by the ratio / tg„ 
where tc is the cooling timescale. The standard cooling 
law due to the modified Urea process is [22J, 

tc{T) w [(T/IO'^K)-^- (T;/10'*K)-^) yr (4) 

where T,- is the initial temperature. Setting T, = 2 x 
10^° K and T = Ten w (5 - 9) x 10*^ K (as suggested by 
the cooling of the neutron star in Cassiopeia A [26 1) we 
obtain f c ~ 1 ~ 100 yr. Thus, the most favourable sce- 
nario where fgw ^ fc is clearly relevant, in which case 
the mode has ample time to grow and eventually satu- 
rate by non-lrnear physics, while emitting gravitational 
waves. 

It is conceivable that newborn neutron stars spinning 
near (1^ are also endowed with strong magnetic fields. 
In this case, the electromagnetic spin-down could op- 
erate on a timescale comparable to the /-mode growth. 
This leads to a timescale roughly given by O /O 
10^ {CiK/Ci)^{{B/10^^G)y^ s. Thus, a magnetar-strong 
magnetic field could effectively choke the instability. 

Concluding discussion — The results presented in this 
Letter add more strength to the possibility that CFS- 
unstable /-modes in fast spinning newborn neutron 
stars could be a promising source of gravitational 
waves. We should, of course, emphasize the obvious 
caveat that the required near-Kepler birth spin rates are 
only expected in a small fraction of newborn neutron 
stars. Our results should serve as the stepping stone to- 
wards a revisit of the gravitational wave observability 
of realistic, rapidly spinning neutron stars. We conclude 
with a discussion of some key issues that should be ad- 
dressed by future work. 

Our analysis made use of a polytropic stellar model as 
an approximation to realistic neutron stars. Future cal- 
culations should explore a wide range of realistic equa- 
tions of state. Although we do not anticipate any qual- 
itative change in our conclusions, some moderate vari- 
ation in the mode's frequency /■ and growth timescale 
fgw is to be expected among different equations of state. 

Another improvement on our work would be to aban- 
don the Cowling approximation and study oscillations 
of rapidly rotating neutron stars in full General Relativ- 
ity. Recent work along this line suggests that the critical 
rotation Cic for the onset of the CPS instability is pushed 
to lower values ||14| , a clear indication of a stronger /- 
mode instability, particularly for the lower-order multi- 
poles. 

The instability could be also enhanced if the stellar ro- 
tation departs from that of a rigid body. Such differen- 
tially rotating stars are known to have a larger instability 



parameter space as compared to the uniformly rotating 
models considered here |9 17J. A key related issue is 
to determine whether differential rotation survives long 
enough as to allow for a significant mode growth. 

It would be also desirable to gain a better understand- 
ing of /-mode damping mechanisms. At the level of lin- 
ear perturbations, it would be worth exploring how the 
instability window is modified if exotic matter (decon- 
fined quarks, hyperons) is present in the neutron star 
core. Although the recent discovery of a 2M0-neutron 
star [27] leaves little room for a hyperonic composition, 
the presence of strange quark matter in the color-flavor- 
locked state remains an open (and exciting!) possibility 



At the non-linear level, it is crucial to set an upper 
limit on the maximum amplitude of the instability. This 
would have obvious implications for the gravitational 
wave observability because the mode could saturate as 
a result of non-linear couplings with other modes. In 
this respect, a recent non-linear study is quite optimistic, 
suggesting observable signals from as far as 10 Mpc 
for Advanced LIGO and VIRGO |29l. Another differ- 
ent non-linear effect that could cause mode saturation is 
non-linear bulk viscosity in the so-called suprathermal 
limit |301. 

Finally, future modeling should account for the pres- 
ence of a crust and its impact on the /-mode instability. 
This could include the damping due to a viscous bound- 
ary layer at the crust-core interface, and even the direct 
breaking of the crust by a large amplitude /-mode with 
the subsequent dissipation of energy at the fracture sites 
ll6ll3Tl. We hope to address some of these exciting issues 
in the near future. 
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